PARAMETER ESTIMATIONS FOR SPDES WITH MULTIPLICATIVE 

FRACTIONAL NOISE 
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Abstract. We study parameter estimation problem for diagonalizable parabolic stochastic partial 
differential equations driven by a multiplicative fractional noise with any Hurst parameter H £ 
(0, 1). Two classes of estimates are investigated: traditional maximum likelihood type estimates, 
and a new class called closed-form exact estimates. Finally several examples are discussed, including 
statistical inference for stochastic heat equation driven by a fractional Brownian motion. 
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1. Introduction 

Parameter estimation problem for stochastic partial differential equation has been of great inter- 
est in the past decade, and besides being a challenging theoretical problem, it finds its roots and 
motivations from various applied problems: fluid dynamics [121 [30], biology [HI [TO], finance [HO [8], 
meteorology [5] etc. At general level the problem is to find or estimate the model parameter -d 
(could be a vector) based on observations of the underlying process u$ which is assumed to be a 
solution of a stochastic evolution equation in finite or infinite dimensional space. We will follow 
traditional continuous time approach and assume that the solution u$(t) is observed continuously 
in time t G [0, T]. From statistical point of view, we suppose that there exists a family of probability 
measures P$ that depends on parameter ■& £ C W 1 , and each P$ is the distribution of a random 
element. Assuming that a realization of one random element corresponds to a particular value i?0) 
the goal is to estimate this parameter from given observations. One approach is to select parameter 
i? that most likely produces the observations. This method assumes that the problem is regular 
or absolutely continuous, which means that there exists a reference probability measure Q such 
that all measures P$, # G O, are absolutely continuous with respect to Q. Then Radon-Nikodym 
derivative dP^/dQ, also called the likelihood ratio, exists, and the Maximum Likelihood Estimator 
(MLE) $ of the parameter of interest is computed by maximizing the likelihood ratio with respect 
to i9. Usually $ ^ $ and the problem is to study the convergence of MLE to the true parameter 
as more information arrives (for example as time passes or by decreasing the amplitude of noise). 
If the measures P# are singular for different parameters then the model is called singular, and 
usually the parameter can be found exactly, at least theoretically. While all regular models are to 
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some extend the same, each singular model requires individual approach. For example, estimat- 
ing the drift coefficient for finite-dimensional stochastic differential equations is typically a regular 
problem, and the parameter can be estimated by means of MLEs, while estimating the diffusion 
(volatility) coefficient is a singular problem and one can find the diffusion coefficient exactly through 
quadratic variation of the underlying process. For some finite-dimensional systems, estimating the 
"drift coefficient" is also a singular problem, and as shown in Khasminskii et al [18] the estimators 
have nothing to do with MLEs. Generally speaking statistical inference for finite-dimensional dif- 
fusions has been studied widely, and there are established necessary and sufficient conditions for 
absolute continuity of corresponding measures (see, for example [23], [ST] and references therein). 
Some of these results have been extended to infinite dimensional systems in particular to parabolic 
Stochastic Partial Differential Equations (SPDE). It turns out that in many cases the estimation 
of drift coefficient for SPDEs is a singular problem, and as general theory suggests one can find the 
parameter "exactly". One of the first fundamental result in this area that explorers this singular- 
ity is due to Huebner, Rozovskii, and Khasminskii |17] , The idea is to approximate the original 
singular problem by a sequence of regular problems for which MLEs exist. The approximation 
was done by considering Galerkin-type of projections of the solution on a finite-dimensional space 
where the estimation problem becomes regular, and it was proved that as dimension of the projec- 
tion increases the corresponding MLE will converge to the true parameter. In [15 ^ 116 1 f25l 126] , the 
problem was extended to a general parabolic SPDE driven by additive noise and the convergence 
of the estimators was given in terms of the order of the corresponding differential operators. For 
recent developments and other types of inference problems in SPDEs see a survey paper by Lotot- 
sky [24] and references therein. Statistical inference for SPDEs driven by multiplicative noise is a 
more challenging problem. First and only attempt to study equations with multiplicative noise is 
given in [6], by considering Wiener (not fractional) type noise without spatial correlation structure. 
Besides MLE type estimates, a completely new class of exact estimates were found due essentially 
to the very singular nature of the problem. 

The aim of this note is to study parameter estimation problem for stochastic parabolic equations 
driven by a multiplicative fractional noise with following dynamics 



where .4,0)*4i an d M. are some known linear operators, W is a fractional Brownian motion with 
a Hurst parameter H 6 (0, 1), and 9 is a real parameter belonging to a bounded set Cl. For 
now, assume that the stochastic integral with respect to fractional Brownian Motion W H is well- 
defined, while the exact meaning will be specified in Section [2.11 The main goal is to estimate the 
parameter 9 based on the observations of the underlying process u(t), t £ [0,T]. Similar problem 
for SPDEs driven by additive space-time fractional noise was investigated in [3 I3T1 [27] . Estimation 
of drift coefficient for finite-dimensional fractional Ornstein-Uhlenbeck and similar processes has 
been investigated by Tudor and Viens [32] for H G (0, 1), Kleptsyna and Le Breton [19] for H 6 
[1/2,1), by developing Girsanov type theorems and finding MLEs. Berzin and Leon [3j estimate 
simultaneously both drift and diffusion coefficients. Least square estimates for drift coefficients 
were established by Hu and Nualart [13], and MLE type estimates for discretely observed process 
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by Hu, Weilin and Weiguo [H]. For a general theory, including Girsanov Theorem and some 
results on statistical inference, for finite dimensional diffusions driven by fractional noise see also 
the monograph by Mishura [28J. 

In this paper we continue to explore the impact of the noise in infinite-dimensional evolution 
equations and its implications on statistical inference. Besides its theoretical roots, this problem is 
also motivated by increasing demand in modeling various phenomena by SPDEs driven by fractional 
noise [5j [11]. We assume that the solution of (P) is observed at every t £ [0,T], and hence each 
Fourier coefficient Uk{t) = (uk(t),hk)H is observable for every t € [0, T], where H is a Hilbert 
space in which the solution leaves and h^, h > 1, is a CONS in H. All results are stated in terms 
of Fourier coefficients In the first part of Section [2] we set up the problem and establish the 
existence and uniqueness of the solution of the corresponding SPDE. In Subsection 12. 21 we introduce 
the main notations and find the MLE for fractional Geometrical Brownian Motion (which is not 
covered explicitly in any other sources, at our best knowledge). In Section [3] we study the estimates 
of drift coefficient 9 of equation (pQ) based on MLE of the corresponding Fourier coefficients. We 
establish sufficient conditions on operators A\ , A± and M. , that guarantee efficiency and asymptotic 
normality of the estimates and some of their versions. Section [5] is dedicated to investigation of a 
new type of estimates called closed-form exact estimates, similar to those studied in [BJ. We show 
that 9 can be found exactly by knowing just several (usually two) Fourier coefficients. Moreover, 
by the same technics we found an exact estimate of the Hurst parameter H too, in both regimes, 
9 known and unknown. Of course there are many other methods of finding the Hurst parameter, 
but it is out of scope of this publication to apply them to our equation. Some of the results follow 
from simple algebraic evaluations, but the very existence of such estimates is amazing and gives a 
better understanding of the nature of the problem's singularity. Also, we want to mention that, 
despite of memory property of the fractional Brownian Motion which is spilled over the solution 
too, the exact estimates are based only on observations at time zero and some future time T. In 
contrast, the MLEs require observation of the whole trajectory u(t), t € [0, T\. We conclude the 
paper with two examples which are of interest along: stochastic heat equation with parameter 9 
next to Laplace operator, and a general second order parabolic SPDE with 9 next to a lower order 
operator. 

While we assume that data is sampled continuously in time, in practice usually this is not the 
case. For the MLEs derived in Section [3] the problem is reduced to approximate some integrals of 
a deterministic function with respect to the solution u and eventually to the fractional Brownian 
motion. However, the Exact Estimates from Section [J] depend only on the values of the solution 
at initial time t = and some future time t = T, and thus do not depend on how the solution is 
observed in time. 

2. Preliminary results 

2.1. The equation and existence of the solution. Let H be a separable Hilbert space with the 
inner product (•, -)o and the corresponding norm || • ||o- Let A be a densely-defined linear operator 
on H with the following property: there exists a positive number c such that ||Au||o > c||u||o for 
every u from the domain of A. Then the operator powers A 7 , 7 € R, are well defined and generate 

the spaces H 7 : for 7 > 0, H 7 is the domain of A 7 ; H° = H; for 7 < 0, H 7 is the completion of H 
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with respect to the norm || • |L := ||A • ||o (see for instance Krein at al. [2D]). By construction, the 
collection of spaces {H 7 , 7 G R} has the following properties: 

• A 7 (H r ) = H r ~ 7 for every y,r G R; 

• For 71 < 72 the space H 72 is densely and continuously embedded into H 71 : H 72 C H 71 and 
there exists a positive number c\ 2 such that ||u|| 7l < ci2||ii|L a for all u G H 72 ; 

• for every 7 G R and m > 0, the space H 7_m is the dual of H 7+m relative to the inner 
product in H 7 , with duality (•, -) 7im given by 

(ui,u 2 ) 1>m = (A 7 - m ui,A 7+m n 2 ) , where m G H 7 " m , u 2 G H 7+m . 
Let (Q, J 7 , {J-t},P) be a stochastic basis with usual assumptions. 

Definition 1. A fractional Brownian motion with a Hurst parameter H G (0, 1) is a Gaussian 
process with zero mean and covariance 

EW H (t)W H (s) = ^(t 2H + s 2H -\t- s\ 2H ), t,s>0. 

Consider the following evolution equation 

du(t) = [(A + 9Ax)u(t) + f(t)]dt + {Mu{t) + g{t))dW H (t), < t < T, 
u(0) = u , 



(2) 



where Aq,A\, A4 are linear operators in H, / and are adapted H- valued processes, uq G H, 
W H is a fractional Brownian Motion with Hurst parameter H G (0, 1), and 9 is a scalar parameter 
bellowing to an open set G C R. 

Definition 2. Equation ([2J is called diagonalizable if the operators Ao, A\ and M have point 
spectrum and a common system of eigenfunctions {hj, j > 1}. 

Denote by pk, Vk, and /ik the eigenvalues of the operators ^4o 5 Ai, and Ai: 

(3) Aoh k = p k h k , Aih k = v k h k , Mh k = p k h k , k > 1, 

and also denote by a k (9) := p k + 0u k , k > 1, the eigenvalues of operator ^lo + 9A\. Without loss 
of generality we assume that the operator A has the same eigenfunctions as operators Aq, Ai, M.: 
Ah k = X k h k , k > 1. 

Definition 3. The equation (|2|) is called parabolic in the triple (H 7+m , H 7 , H 7_m ), for some positive 
m and real 7, if there exists positive real numbers 5, C\ and a real number C 2 such that, for all 
k > 1 and all G G, 

(4) \ k 2m \Pk + 9u k \ <C i; 

(5) 2(p fc + ^ fc )+/i2+^ m <C 2 . 

This definition is equivalent to the classical definition of parabolic equations, but written in terms 
of eigenvalues of corresponding operators. 

Theorem 4. Assume that equation ([2]) is diagonalizable and parabolic in the triple (H 7+m , H 7 , H 7_m ), 
the initial conditions uq is deterministic and belongs to H 7 , the process f = f(t) is Tt-adapted with 



values inWi~ m , andE \\f(t)\\ 2 _ m dt < oo, the process g = g(t) is Tt-adapted with values in H 7 
and E J Q T ||<?(i)||^ dt < oo. Then the process u defined by 

(6) u(t) = J2 u k(t)hk, 

k>i 

where 

(7) u k (t) = u fc (0)exp ([a k (6) + f k (t)]t - \[^k + 9k(t)] 2 t 2H + \pt k + g k {t)]W H (t) 

/(*) = E g{t) = Y,9k{t)h k , 

k>l k>l 

is an H 7 -valued stochastic process. 

Proof. Since W H (t) is a Gaussian random variable with zero mean and variance t 2H , we have 
K\u k (t)\ 2 = ul(0)exp {2{a k {6) + f k (t))t + (fj, k + g k {t)) 2 t 2H ) . 

Hence, 

E|| U (t)B = ^A^K(i)| 2 < C^exp(2a fc (^ + ^ 21f ) . 

fe>l k>l 

By parabolicity condition ([5]), the last series converges uniformly in i, and the theorem follows. □ 

The functions formally represent the Fourier coefficients of the solution of equation (j2]) with 
respect to the basis {h k } k >i and the uniqueness of u follows. Since the equation is diagonalizable, 
naturally we conclude that formally u k has the following dynamics 

(8) du k (t) = (9u k + p k )u k {t)dt + fi k u k {t)dW H (t), k > l,t > 0. 

Specifying the stochastic integration in ([2]) is equivalent to specifying in what sense we understand 
the integration with respect to fractional Brownian Motion for the Fourier coefficients ([8]). Con- 
sequently, since the equation has constant coefficients, specifying the solution of (|8|) is the same 
as to stipulate the sense of stochastic integration in ([8]). If the integration is understood in Wick 
sense then u k , k > 1, defined in ((JJ) is the unique solution of equation (jHJ) for all H E (0, 1) (see 
for instance [4], Theorem 6.3.1). All results stated here are easily transferable to any other form of 
integration, by caring out the relationship between different form of integration and consequently 
adjusting the form of the solution of equation ([8]) (for comparison of various form of integration 
with respect to fBM see 0], Chapter 6). Our choice was just to have a unified theory and same 
formulas for all H € (0, 1). 

Definition 5. The process u constructed in Theorem |4] is called the solution of equation ([2]). 

It should be mentioned that the above result, with some obvious adjustments, also holds true 
for diagonalizable equations driven by several independent fractional Brownian Motions, even with 
different Hurst parameters. 



2.2. Parameter estimation for geometrical fractional Brownian motion. In this section 
we will present some auxiliary results about parameter estimation for one dimensional diffusion 
processes driven by multiplicative fractional noise. For similar results for equations with additive 
noise see for instance Kleptsyna and Le Breton |19j . Tudor and Viens [32] . or Mishura [28], Chap- 
ter 6. The results essentially follow from Girsanov type theorem for diffusions driven by fractional 
Brownian motion. 

Let r and B denote the Euler Gamma- functions. Following Mishura [28] we introduce the 
following notations 

= / T(3 - 2H) 

H \2HT(^-HfT{\ + H) 

(9) l H (t, s) = C H s l 2- H (t - s)^ H I 0<s<t , 

Mf:= [ t l H (t,s)dW s H , , 
Jo 

where H 6 (0, 1), and the integration with respect to fractional Brownian Motion is understood in 
Wiener sense (for more details see [28], Chapter 1). The process is a martingale, also called 
the fundamental martingale associated with fractional Brownian motion W t H (see for instance [29] 
or [28], Theorem 1.8.1). has quadratic characteristic (M H )t = t 2 ~ 2H , and by Levy theorem, 
there exists a Wiener process {Bt,t > 0} on the same probability space such that 

M t H = (2 - 2H)2 / S 2- H dB s . 
Jo 

Moreover, a{W s H ,0 < s < t) = a(B s ,0 <s<t). 
Let us consider the stochastic process of the form 

X t = X exp (et - ^ 2 t 2H + aW H {t)J , t > 1 , 

which can be called the Geometric Fractional Brownian Motion, and as mentioned in the previous 
subsection it is the unique solution of the stochastic equation 

dX t = eX t dt + aX t dW t H , Xq = xo, t G [0, T). 

Let Y t := lnX t /X = Ot - + aW t H , and consider the process Y t := f*l H (t,s)dY s . Note 

that observing one path of the process {Y s , < s < t} implies that the one path of process 
{Y s , < s < t} is observable too. By Q we have 

(10) Y t = aM t H + eb x t 2 ~ 2H - a 2 Hb 2 t, t > , 

where hi = C H B(3/2 - H, 3/2 - H), b 2 = C H B(l/2 + H, 3/2 - H). 

For a fixed parameter 6 £ 0, let us denote by Pg the distribution of the process Yt and by Po 
the distribution of the process Y® := aM^ 1 = ab^ s l l 2 ~ 2H dB s . The measure P# is absolutely 
continuous with respect to Po and the Radon-Nikodym derivative, or the likelihood ratio, has 
the following form (see for instance [23J, Theorem 7.19 or apply classical Girsanov Theorem for 
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martingales) 

dF e ( f °( 2 ~ 2H)b lS l - 2H - a 2 Hb 2 ~ , 1 f* [9(2 - 2H)b 1 s l ~ 2H - a 2 Hb 2 } 2 

dF {Yt)=eW {-Jo dYs + lJo ^ 

The MLE is obtained by maximizing the log-likelihood ratio with respect to 9. Since 

8 dF e .~* (2-2H)h~ i (2-2H)b\t 2 - 2H (2-2H)b l Hb 2 t 

ae ln dF- { t) ~ ^bV~ * + s% 6g ' 

the MLE for parameter 9 has the form 



b X t 2 ~ 2H ' bit 1 ' 211 ' 

Proposition 1. The estimate 9t, t > 0, is an unbiased estimate for parameter 9q ; lim 9t = 9q with 

t— >-oo 

probability one, i.e. 9t is a strong consistent estimate of 9q; t 1 ~ H (9t — 9q) converges in distribution 
to a Gaussian random variable with zero mean and variance a 2 jb\. 

Proof. Using the definition of the process Yt, we represent the estimate 9t as follows 
(12) 9 t = 9 0+ ° M ? 



ht 2 - 2H ' 
where 9q is the true parameter. 

The unbiasedness and asymptotic normality follows immediately from (|12p and the fact that 
is a Gaussian random variables with zero mean and variance t 2 ~ 2H . Since is a square 
integrable martingale with unbounded quadratic characteristic t 2 ~ 2H — > oo, as t — > oo a.s., by 
Law of Large Numbers for Martingales [22], Theorem 2.6.10, j (M H )t — > a.s., and hence 
consistency follows. □ 

Note that, in particular, for H = 1/2 we have b\ = b 2 = 1, and we recover the classical estimate 
for the drift coefficient of geometric Brownian Motion 

e t - J + -- — - — + __e + ^, t >o, 

and its corresponding asymptotic behavior. 

3. Maximum Likelihood Estimator for SPDEs 
Consider the diagonalizable equation 

(13) du(t) = (A + Ai)u(t)dt + Mu(t)dW H (t), 

with solution u(t) = Y2k>i u k(t)hk given by ([7]). As mentioned in Introduction, if u is observ- 
able, then all its Fourier coefficients Uk can be computed. Thus, we assume that the processes 
u\(t), . . . ,UN(t) can be observed for all t € [0, T] and the problem is to estimate the parameter 9 
based on this observations. Also, we assume that the Hurst parameter H G (0, 1) is known for now. 

By Definition [5] of the solution of equation (|13p the Fourier coefficients Uk, k £ N, have the 
following dynamics 

(14) dukit) = a k (9)u k (t)dt + n k u k (t)dW H (t), t £ [0, T], 
where otk{9) = pu + 9vk, k G N. 



For every non-zero ttfc(O), k € N, denote by ffc(i) = ln(nfc(t)/ufc(0)), and Sfc(t) = l(t, s)dvk(s), 
where /(•,•) is denned in ([9]). By results of Section 12.21 it follows that there exists a Maximum 
Likelihood Estimate for ak(6) and it has the form 

w> 55) = ^+^, '>>■ 

Since a.k(9) is a strictly monotone function in 0, by invariant principle of MLE under invertible 
transformations, we can find an MLE for the parameter 9 

(16) ^ = ^&r + 4^"-' *^Me[o,H. 

Using the definition of the process Vk, the estimate Ok,T can be represented as follows 
(in ft -ft i ^ m t 

and by similar arguments to the proof of Proposition [H we have the following result. 

Theorem 6. Assume that equation f)13|) is diagonalizable and parabolic in the triple (H 7+m , H 7 , H 7 ~ 
for some -y € R, m > and no € H 7 . Then, 

(1) For every k > 1 and T > 0, is an unbiased estimator of 9q. 

(2) For every /ixed > 1, as T —?■ oo, the estimator 9k,T converges to 9q with probability one 
and T 1 ~ H (9k i T — @o) converges in distribution to a Gaussian random variable with zero 
mean and variance ^\/b\u 2 . 



(3) If, in addition, 



(18) lim 

k— >oo 



0, 



i/ien for every fixed T > 0, lim 0^ = #o with probability one, and \uk/ Hk\{9k,T — Oo) con- 
verges in distribution to a Gaussian random variable with zero mean and variance T 2H ~ 2 /b\. 

Remark 7. The parabolicity conditions (U])-© and MLE consistency condition (|18|) in general are 
not connected. In terms of operator's order, parabolicity states that the order of operator M from 
the diffusion term is smaller than half of the order of the operators Aq and A\ from deterministic 
part. Condition (|18p . that guarantees the consistency of MLE as number of Fourier coefficients 
increases, assumes that the order of operator M. from the diffusion part does not exceed the order 
of the operator A\ from deterministic part that contains the parameter of interest 9. 

By Theorem [6] it follows that the consistency and asymptotic normality of the estimates 6k,T can 
be achieved in two ways: by increasing time T or by increasing the number of Fourier coefficients 
k. In both cases the quality of the estimate is improved by decreasing its variance. 

It is interesting to note that Var \ 9k,T — = h L 'k^ 2H ' 2 /^i h 'k a ^ so depends on Hurst parameter 
H. For H > 1/2 the constant l/b± is close to one, and increases as function of H for H £ (0, 1/2). 
The function t 2H ~ 2 increases in H for any t > 1. The constants u^. and Uk, k > 1, do not depend 
on H. Overall, T 2H ~ 2 /b 2 increases in H for any t > 1 and thus quality of the estimates is higher 
for smaller H. 



8 



As mentioned before, due to the fact that the probability measures generated by the solution u 
of the original SPDE are singular, it is possible to estimate 9 exactly on any finite interval of time 
[0,T]. A natural question is wether we can improve the quality of the estimates by considering 
several Fourier coefficients Uk(t). The answer is that by statistical methods used above this is not 
possible. First, note that the measures associated to any two or more processes are singular, 
and thus MLE does not exist for such vector- valued functions. In other words, by considering two 
or more Fourier coefficients Uk, we get a singular model, a fact that will be explored in the next 
section. Also, since each process Uk is driven by the same noise, each individual Fourier coefficient Uk 
contains the same amount of information: the sigma-algebra generated by Ufc(t), t E [0,T] coincides 
with the sigma-algebra generated by W H (t), t £ [0, T]. However, the speed of convergence of the 
sequence 9kT can be improved by using accelerating convergence technics from numerical analysis. 
Two methods have been discussed into details in [6]: the weighted average method and Aitken's 
A 2 method. For sake of completeness, we will state here the corresponding results applied to the 
sequence {0k,T}k>i- 

Weighted averaging. Suppose that /3k, k > 1, is a sequence of non- negative numbers such that 
Sfe>i Pk = +oo, and consider the weighted averaging estimator 

N N 

(19) 

"(N,T) — ^Pkhp/^Pk N>1,T>0. 
k=l k=l 

Then (a) 0(nt) IS an unbiased estimator of On f° r every N > 1 and T > 0; (b) lim 9( NT -. = $o 
a.e. for every N > 1 (consistency in T); (c) if in addition the consistency condition (|18p is fulfilled, 
then lim 9(nt) = with probability one for every T > (consistency in N). 

JV-s>oo ^ ' ' 



Aitken's A 2 method. Define the following sequence of estimates 



(20) 



k+2,T + 29k+l,T — &k,T 



One can show that the new sequence 6^ converges to the true parameter #o with probability one. 
Moreover, if fik/^k ~ ak~ s for some a, 5 > 0, then 

M9k - n ) 2 1 



and if \ikl v k = {~^) k /k, then 



~ ji), c>0. 

In both cases, the new sequence 6^ converges faster than to #o- 

The proofs of the above results follows from Theorem [6] and some direct computations, which 
will be omitted here. 
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4. Exact Estimates 



In regular models the unknown parameter can be found only approximatively, and the consistency 
is gained either in large sample or small noise regime. For singular models the parameter can be 
found exactly. For example, if all Fourier coefficients of the solution u of equation ([2]) are known, 
according to the results from previous sections, one can find the value of 9q exactly, on any interval of 
time [0, T]. The possibility to evaluate 9q exactly is based on singularity of the measures generated 
by u e for different values of 9. However, while theoretically it is possible to estimate the true 
parameter exactly, in practice we (or computer) can perform only a finite number of operations. 
Recall that the measures associated to an individual Fourier coefficient v? k are regular, while a 
vector consisting of any two or more Fourier coefficients will produce measures that are singular. 
In this section we will explore this singularity, and show that in fact the true parameter can be 
estimated exactly from a finite number of Fourier coefficients. Moreover, the described method 
allow to find both parameters 9 and H, either individually or simultaneously. 

Following [6] we say that an estimator is closed-form exact if it produces the exact value of the 
parameter of interest after finite number of additions, substraction, multiplications, and divisions 
on the elementary functions of the observations. 

Closed- form exact estimates exists for the model ([2]) if we assume that observations are u k (t), k > 
1, t G [0, T\. For every non-zero Fourier coefficient u k of the form (fl4j) . set Vk(t) = lnu k (t)/u k (0), t € 
[0,T]. Then 

(21) v k (t) = (p k + 9v k )t - X -p\t 2H + p k W H (t) . 

Case 1. 6 unknown, H known. Assume that v k p m 7^ v mPk f° r some k, m € N. Then, taking (|21|) 
for these k and m, by direct arithmetic evaluations, one gets the exact estimate of the parameter 9 

, 00 , „ PmVk ~ PkVrn + {PmPk ~ PkPm)t + \{p 2 k Prn ~ P m Pk)t 2H 

(22) a = , 

t[V k Prn — VmPk) 

for any t > and k, m € N for which ffc// m 7^ Vrn^k- 

Note that if ^ = /i m then the above exact estimate does not depend on H, and 9 can be 
evaluated even if H is unknown. This is the case, for example, if M is the identity operator (see 
Example 1 below). 

Case 2. H unknown, 9 known. Assume now that the parameter of interest is the Hurst 
parameter H and assume that 9 is known. By the same arguments as above, one can solve for 
H the system of two equations generated by (|2ip for some k and m, and get the following exact 
estimate for H 

(Pk + 9v k )t[i m - (p m + 9v m )t\i. k - v k p m + V 

mpk 

pm ) 

for any t > 0, k ^ m, and under assumption that the expression under logarithm is positive and 
finite. 

Case 3. Both 9 and H unknown. Denote by a fc>m := {v k p m -v m Pk)t, Pk,m ■= ^/^(p m Pk-p 2 k Pm) 

and4,m VkPm-VmPk- PkPmt- PmPkt- Assume that for some m,k,i,j G N, a k , m (3 itj ^a k , m (3 it j. 
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(23) 



H 



21nt 



In 



Then the following exact estimate for holds true 



(24) 



.in 

If in addition 6k,m<^i,j 3i,j a k,m, then there exists an exact estimate for Hurst parameter H given 
by 

(25) H=^\og 



1 , fok^mQij $i,jCtk,m 
" Pk,m^i,j Pi,j(^k,m 

Note that for this case, generally speaking, it is sufficient to know only three Fourier coefficients, 
i.e. some of the indices k,m,i,j can coincide. 

Remark 8. 

(a) Applying the above idea, closed-form exact estimates can be obtained for equations driven 
by several fractional Brownian motions, even with different Hurst parameters. If we assume 
that the noise is driven by n fBMs, and that one of the parameters 9 or H is known, 
then by considering n + 1 Fourier coefficients we can eliminate all noises and get a closed- 
form estimate as a solution, under some non-degeneracy assumptions. Respectively, if 
both parameters are unknown, then one can estimate them by considering n + 2 Fourier 
coefficients. 

(b) Note that the construction of the exact estimates assumed only the existence of the solution 
and did not impose any additional assumptions on the order of the operators Ao, Ai, M, in 
contrast to MLE estimates where the consistency holds only under additional assumptions 
on order of corresponding operators. 

(c) The MLE 0j.t depend on the whole trajectory of the Fourier coefficient Uf;(t), t £ [0,T]. 
All exact estimates depend only on initial and terminal value of life's. 



5. Examples 

We conclude the paper with two practical examples where we explore some of the estimates 
proposed above. 

Example 1. Stochastic heat equation. Let 9 be a positive number, and consider the following 
equation 

(26) du(t, x) = 9u xx (t, x)dt + u(t, x)dW H (t), t > 0, x 6 (0, 1), 

with zero boundary conditions and some nonzero initial value u(0) € £2(0, 1). In this case the oper- 
ator Ai is the Laplace operator on (0, 1) with zero boundary conditions that has the eigenfunctions 
hk{x) = y2/7r sin(/cx), k > 0, and eigenvalues = —k 2 , = 0, = 1, k > 0. Assume that 
u(t,x) is known for x G [0, 1] and t € [0, T], hence Ufc(i) := J hk{x)u(t,x)dx, k 6 N, is observable. 
Denote by Vk(t) := log(uk(t) / Uk(0)) for every k G N, and u^(0) ^ 0. By Theorem[6l the MLE for 
9 has the form 

Iol(T,s)dv k (s) Hb 2 

k ~ k 2 hT 2 - 2H k 2 b x T^- 2H ' ^ 
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The exact estimates (|22|) for 9 are given by 



1 ^u k (T)u m (0) 



T(m 2 -k 2 ) u m (T)u k (0)' 

for any k ^ m and T > 0. Note that the exact estimates do not depend on H. However there are 
no exact- type estimates for H. 

Example 2. Assume that G is a bounded domain in M. d , and let A be the Laplace operator on G 
with zero boundary conditions. Then A has only point spectrum with countable many eigenvalues, 
call them o~ k , k £ N. Moreover, the set of corresponding eigenvalues forms an orthonormal basis 
in L2(G); the eigenvalues can be arranged so that < —o\ < — 02 < ■ ■ ■; the eigenvalues have 
the asymptotic a k ~ k 2 l d . In the space H°(G) let us consider the following stochastic evolution 
equation 

du(t) = [Au(t) + 9u{t)]dt + (1 - A) r u(t)dW H {t) , 

with some nonzero initial values in H°(G), and some r € BL According to our notations we have 
the operators Aq = A, A\ = I, M. = (1 — A) r , with corresponding eigenvalues v k = 1, p k = 
°~ki Mfe = (1 + °7c) r - The equation is diagonalizable, and by Theorem [H it has a unique solution in 
the triple (H 1 ,^,^ 1 ) for any r < 1/2. 

The maximum likelihood estimate in this case has the form 
2- v k {t) Hb 2 (l-a k ) 2r 1 

which is an unbiased estimate of the parameter 9. 

(2a) Large time asymptotics. lim 6 k t = 0q a.s. for all k > 1; lim t 1 ~ H (9 k t — 60) = where 

t— >oo ' t—^oo 

£~M(0,(l-a k ) 2 /bj). 

(2b) Consistency in number of spatial Fourier coefficients. Assume that r < 0. Then lim k t = #0 

k— ¥00 ' 

a.s., for every t > 0, and the sequence (1 — o~ k )~ l {9 k)t — 9q) converges in distribution to a Gaussian 
random variable with mean zero and variance t 2H ~ 2 /b 2 . If r € [0,1/2] the solution still exists in 
the space H°(G), while the estimate 9 k> t is not consistent in k. 

(2b) Exact estimates. Let v k (t) = ln(u k (t)/u k (0)). Assume that Hurst parameter H is known. 
Then we have the following exact estimate for 9 

a _ (1 - o- m ) r v k - (1 - o k ) r v m a m (l - a k ) r - a k (l - a m ) r 



t((i - a m y - (1 - a k y) (1 - a m y - (i - 

t 2H ~ l (1 - a k ) 2r (l - a m y - (1 - a m ) 2r {l - a k ) r 
2 (1 - a m y - (1 - a^r 

for any k ^ m and £ > 0. 

If 9 is known, then the Hurst parameter H can be found by 

1 , [(a k + 9){l - a rn ) r - (a m + 9){l - a k ) r ]t - v k (l - a k y + v m (l - a Ti 



H = -log 



2(1 -a k ) 2r (l-a m y -(l-a m ) 2 -(l-a k y 
for any k ^ m, t > 0. 

Finally one can write the exact estimates (|24p and (|25p for the case when both parameters 9 and 

H are unknown. Note that the exact estimates exists for all r as long as the solution exists (maybe 

in a larger space) and the Fourier coefficients u k (t) are computable. 
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